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^ Abstract 

I""! We have investigated the magnetoelectric and magnetodielectric response in FeV04, which ex- 

t/3 hibits a change in magnetic structure coincident with ferroelectric ordering at T/v2~15 K. Using 

\~i symmetry considerations, we construct a model for the possible magnetoelectric coupling in this 

. system, and present a discussion of the allowed spin structures in FeV04. Based on this model, 

rj in which the spontaneous polarization is caused by a trilinear spin-phonon interaction, we experi- 

1-^ mentally explore the magnetoelectric coupling in FeV04 thin films through measurements of the 

Q electric field induced shift of the multiferroic phase transition temperature, which exhibits an in- 

O 
'~~' crease of 0.25 K in an applied field of 4 MV/m. The strong spin-charge coupling in FeV04 is also 

►v^ reflected in the significant magnetodielectric shift, which is present in the paramagnetic phase due 

O 

■^ to a quartic spin-phonon interaction and shows a marked enhancement with the onset of magnetic 

cn 

■^ order which we attribute to the trilinear spin-phonon interaction. We observe a clear magnetic 

^-^ field induced dielectric anomaly at lower temperatures, distinct from the sharp peak associated 



o 



with the multiferroic transition, which we tentatively assign to a spin reorientation cross-over. We 
also present a magnetoelectric phase diagram for FeV04. 
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C^ PACS numbers: 75.85.+t, 77.55.Nv, 75.25.-j 



I. INTRODUCTION 

Magnetoelectric multiferroics, insulating magnets exhibiting simultaneous magnetic and 
ferroelectric order, are widely investigated in large part due to their potential applications 
for developing novel devices, including magnetic sensors and multistate memory, among 
others pl2] . The cross-control of these distinct order parameters, such as adjusting the 
magnetization using an applied electric field or vice-versa, is expected to provide an extra 
degree of freedom in developing new types of spin-charge coupled devices, such as voltage 
switchable magnetic memories [3] . Additionally, there are a number of fundamental materi- 
als questions surrounding the development of multiferroic order. Magnetic and ferroelectric 
order are generally contraindicated in the same phase, as ferromagnetism in transition metal 
systems typically requires partially filled d-orbitals while ferroelectric distortions are pro- 
moted in a d° electronic configuration^. Despite this apparent restriction, a rather large 
number of single phase systems have been identified as magnetoelectric multiferroics [5H7] . 
A number of microscopic mechanisms have been proposed for the development of multi- 
ferroic order, including, a magnetic Jahn- Teller distortion[S] for TbMn205,[n] bond and 
site ordering having distinct centers of inversion symmetry, [lOj a microscopic mechanism 
leading to a spin-current interaction [11], the Dzyloshinskii-Moriya interaction, |[12] a gen- 
eral anisotropic exchange striction,^3] a spin-phonon interaction, [14j and a strain induced 
ferroelectricity. [T3] 

Phenomenologically, magnetically-induced ferroelectric order developing in systems hav- 
ing multiple magnetic phases can be understood by considering a trilinear term in the 
magnetoelectric free energy, Fme, coupling the electric polarization with two distinct order 
parameters ai and a2 which together break inversion symmetry so that Fme ocPcti cq plTllTGl - 
[18] . Since the free energy must transform as a scalar, there are strong symmetry restrictions 
on the allowed representations for ai and a2] in particular, the product cri(q)cr2(q)* must 
be antisymmetric under spatial inversion. This trilinear coupling also predicts electric field 
control of a magnetic order parameter p^lTO] . A general discussion of the symmetry of the 
magnetoelectric coupling in multiferroics is considered for the specific case of FeV04 in the 
following section. Investigations on multiferroic Ni3V208 thin films, in which such a trilinear 
coupling is believed to be responsible for the multiferroic order [Tf7] . have established that 
the multiferroic transition temperature can be varied through the application of either (or 



both) magnetic and electric fields [19] , confirming the strong coupling between magnetic and 
dielectric degrees of freedom. Higher order magnetoelectric coupling terms quadratic in both 
magnetic and ferroelectric terms will give rise to magnetization induced shift in the dielectric 
response[20]. Such coupling has been investigated both theoretically and experimentally in 
a range of materials including Mn304|14j, CoCr204[2l], BaMnF4 [22H25] , and SeCuOs and 
TeCu03[26|. Because this magnetodielectric coupling is also expected to depend strongly 
on the symmetry of the magnetic phase, it has been suggested that changes in this coupling 
may be used to probe changes in the ordered spin structure [H]. 

Triclinic iron vanadate, FeV04, has recently been identified as a multiferroic system 
having the PI space group [27H29] . Magnetic, thermodynamic, and neutron diffraction 
studies on FeV04 single crystal and ceramic samples have shown that FeV04 transitions 
from a paramagnetic phase into a collinear incommensurate (CI) phase at T/vi=22 K and 
then into non-coUinear incommensurate (NCI) phase at T/v2=15 K. [27H30] Ferroelectric order 
in FeV04 develops in this non-collinear spiral magnetic phase. The onset of ferroelectric order 
with the development of a second magnetic phase suggests that a symmetry-based approach 
may be useful in exploring the multiferroic properties in this system. We present a full 
Landau theory for this system, specifically considering the allowed magnetoelectric coupling 
terms. Our result is that a nonzero induced spontaneous polarization P requires having a 
magnetic spiral [31] described by two order parameters which are out of phase with respect 
to one another. [T|7]ll(j||18] In this low symmetry structure there are no restriction on the 
orientation of P based on symmetry arguments, unlike the majority of similar magnetically- 
induced multiferroics [T]l71IT6]l3T] . 

This paper is organized as follows. In Sec. II we present a symmetry analysis of FeV04 
based on Landau theory. Here we analyze the symmetry of the magnetoelectric interaction. 
In Sec. HI we present the results of a number of experiments designed to probe the structure 
of these magnetoelectric interactions. In Sec. IV we briefiy summarize our results. 

II. LANDAU THEORY 

Motivated by this general discussion of the possibility of magnetically driven ferroelec- 
tric order in FeV04, we now present a Landau theory for FeV04 with some details of the 
construction relegated to the the Appendix. As discussed in detail in Ref. [H], the Fourier 



transform of the spin ordering is proportional to the critical eigenvector of the inverse sus- 
ceptibility matrix at the ordering wave vector q. In the Appendix we analyze the constraint 
of spatial inversion in the PI space group of the paramagnetic phase, with the following 
results. The FeV04 structure consists of six S=5/2 Fe^"*" spins in the unit cell at locations 
Tn- For n = 1, 2, 3, -r^ = r „ = r„+3 and -t„+3 = t„+3 = r„. Then inversion symmetry 
(X) implies that spin Fourier transform obeys 

IS{q,T) = S{q,ry (1) 

where, as defined in the Appendix, S{q, r) is the spatial Fourier transform of the thermally 
averaged spin operator. As explained in the Appendix, this relation implies that the spin 
distribution is inversion-symmetric about some origin. We find that inversion symmetry 
implies that 

[S.{l),Sy{l),S,{l),S.{2),Sy{2)...S.{Q)] = a4xl,ylzl,x;,y;,...4] (2) 

where all the components are complex valued with Xn = x^, yn = 1/^, Zn = z^, are normahzed 
by Z]n=i[|2;nP + Iz/nP + knP] = 1, and the wave vector argument is implicit. The amplitude 
CTniq) is the complex valued magnetic order parameter, which obeys 

As noted in the Appendix, this relation implies that each o"„ is inversion invariant about a 
lattice point (which depends on n) where the order parameter wave has its origin. As the 
temperature is lowered one passes from the paramagnetic phase into a phase with an order 
parameter cri(g) and then, at a lower temperature into a phase where two order parameters 
(Ji{q) and (J2{q) are nonzero, both of which obey Eq. (|3]), but which have different centers 
of inversion symmetry. 

The total magnetoelectric free energy, Fme can be written as 

FuE = Fu + Fe + V , (4) 

where Fm is the purely magnetic free energy, F^ is the dielectric potential which we ap- 
proximate as Fe = (l/2)x^^P^, where xe is the dielectric susceptibility (whose crystalline 
anisotropy is neglected), and to leading order in an, the magnetoelectric coupling term is 
given by: 

2 

V= J2 J2i(^n,'n,jcrn{q)crm{q)* + a*^„^^^an{q)*a^m{q)]P^ , (5) 

n,m=l 7 



where n and m label order parameter modes and 7 labels the Cartesian component of P. 
Terms linear in an are prohibited because they are not time reversal invariant and also can 
not conserve wave vector. The magnetoelectric interaction V has to be inversion invariant 
and the appendix shows that the a coefficients are pure imaginary, so that 

V = z^r^[ai(g)(T2(g)* -(7i(g)V2(g)]PT, = 2^r^|ai(g)CT2(g)|sin(02-0i)A' ' (6) 

7 7 

where o"„(g) = |cr„(g)| exp(z</)„). There is no restriction on the direction of the spontaneous 
polarization, so that all components of P will be nonzero. However, if the magnetic structure 
is a spiral, then the arguments of Mostovoy[3T] might be used to predict the approximate 
direction of P. The result of Eq. (6) is quite analogous to that for Ni3V208[7] or for 
TbMnOalin], in that it requires the presence of two modes cri{q) = exp{i(f)i)\ai{q)\ and 
(j2(g) = exp{i(f)2)\o'2{q)\ which are out of phase with one another: 0i 7^ 02- Then the order 
parameter wavefunctions have different origins and will therefore break inversion symmetry. 
If, as stated in Ref. EHl the eigenvector is not inversion invariant as implied by Eq. (Is]), 
then one would conclude that the magnetic ordering transition is not continuous. However, 
the most likely scenario is that the ordering transitions are continuous and that the spin 
distribution for each an{(j) is inversion symmetric as obtained in this derivation. The acentric 
distribution found in Ref. 128) differs only slightly from being inversion symmetric for reasons 
that are obscure. 132) 



III. MAGNETOELECTRIC INTERACTIONS (EXPERIMENTAL) 

A. Sample Synthesis and Structural Characterization 

Motivated by Eq. (6), which predicts that the magnetic structure defined by cri(g) and 
(T2{(j) is coupled to the electric polarization P, we experimentally investigated the nature 
of the higher order magnetoelectric coupling in FeV04. Bulk single phase polycrystalline 
iron vanadate (FeV04) ceramic samples were prepared using standard solid state reactions. 
Because Eq. 6 predicts that all components of the polarization vector are nonzero, and pre- 
vious measurements on ceramic FeV04 have found clear evidence for multiferroic behaviour 
[27] we focused our study on polycrytalline samples. A stoichiometric ratio of iron oxide 
(Fe203) and vanadium pentaoxide (V2O5) solid solutions were thoroughly mixed and ground 



to produce a homogeneous mixture. This mixture was slowly heated to 600°C for 4 hours in 
air. Intermediate grindings followed by thermal annealing in air were repeated several times 
to complete the solid state reaction and ensure a fully reacted and uniform composition. 
This homogeneous solid solution was finally annealed in air at 800°C for 4 hours, yielding a 
yellowish brown powder identified as a single phase iron vanadate by X-ray diffraction and 
Raman spectroscopy. 

In order to apply large electric fields to FeV04 we also prepared thin film samples. These 
were fabricated from a phase pure stoichiometric iron vanadate target. The FeV04 powder 
used for the sputtering target was prepared by the method described above. Approximately 
30 g of FeV04 powder was mixed with 15 mL of 2 mole percent polyvinyl alcohol as a binder. 
The dried powder was pressed into a circular disc having a diameter of approximately 50 mm 
with a thickness of roughly 3.5 mm followed by air annealing at 600°C for 4 hours to burn off 
the residual organics. A final thermal annealing was done at 800°C for 4 hours to produce the 
dense pellet used for the sputtering target. FeV04 films were deposited at room temperature 
using RF magnetron sputtering onto conducting silicon substrates. The working pressure 
was held at 1.5x10^^ torr, with the atmosphere consisting of a mixture of approximately 
1.5x10^^ torr partial pressure of oxygen as the reactive gas and approximately 1.35x10"^ 
torr partial pressure of argon as the sputtering gas. These as deposited films, prepared over 
a time of 4 hours, were amorphous. After air annealing at 700°C for 4 hours the films were 
indexed as single phase poly crystalline FeV04. 

We investigated the structural, magnetic, and electronic properties of these samples using 
a number of different techniques. We used a Rigaku RU200 powder X-ray diffractometer and 
Horiba Triax Raman spectrometer to study the crystalline structure of these samples. We 
used a Hitachi scanning electron microscope (SEM) to investigate the surface morphology 
of the thin film samples and an associated energy dispersive X-ray (EDX) assembly to 
probe the chemical composition of both samples. We measured the temperature dependent 
magnetization of the powder sample using a Quantum Design MPMS SQUID magnetometer, 
although the very small magnetic anomalies associated with the transitions could not be 
clearly distinguished from background in the thin films samples. We conducted temperature 
and field dependent dielectric and pyrocurrent measurements using the temperature and 
field control provided by a Quantum Design PPMS system used in conjunction with an 
Agilent 4284A LCR meter and a Keithley 6517 electrometer. These measurements were 
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done on a cold pressed pellet of bulk FeV04 with the top and bottom electrodes fashioned 
using silver epoxy and on the FeV04 thin films with room sputtered gold (Au) used as the 
top electrode and the Si substrate serving as the bottom electrode. 

The structure of the ceramic FeV04 sample was practically identical to that previously 
presented for a bulk sample prepared using a different technique |27i|. The X-ray diffraction 
(XRD) pattern for the FeV04thin film is shown in Fig 1(a). These diffraction peaks are 
consistent with the expected XRD pattern for FeV04 [JCPDS no. 38-1372]. The surface 
morphology the thin film sample is shown in Fig 1(b). This SEM micrograph indicates 
that the film consists of grains with various orientations, as well as a number of pinhole 
defects. We calculated the thickness of these thin films to be roughly 200 nm, using the 
cross-sectional SEM micrograph, Fig 1(c). This value is very consistent with estimates from 
well defined interference fringes observed in reflection spectra (not shown). EDX analysis 
of both the bulk and thin film samples show a 1:1 iron to vanadium ratio. We carried out 
room temperature Raman vibrational spectroscopy to further probe the microstructures of 
both bulk and thin films. The identification of Raman active modes and their detailed 
temperature dependent analysis on bulk FeV04 sample is discussed elsewhere [27] . Here we 
plot the room temperature Raman spectrum of both bulk and thin film FeV04in Figl(d). 
We are able to identify all the Raman active modes for thin films, which are observed in 
bulk FeV04 [27] with a small shift in the Raman peaks for the thin films. The Raman peak 
arising from the silicon substrate is indicated by an asterisk. 

B. Temperature Dependent Dielectric Measurements on FeV04 Ceramic 

The temperature dependent magnetization for the bulk FeV04 sample (not shown) was 
practically identical to that measured previously on a different ceramic sample prepared 
using a different technique [27]. In particular, the magnetization showed the usual two 
anomalies associated with the two incommensurate transitions in this system. We plot the 
zero-field dielectric constant and resistivity for bulk FeV04 over a broad range of temper- 
atures in Fig. 2(a). The dielectric constant exhibits a sharp peak near T/v2, arising from 
the development of ferroelectric order in the incommensurate spiral magnetic phase. As 
shown in Fig. 2(a), above 35 K, the dielectric constant for FeV04 shows a gradual decrease 
on cooling, typical of many insulating materials |33j. Below roughly 30 K, the dielectric 
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FIG. 1: (a) 0-29 x-ray diffraction (XRD) pattern of FeV04 thin film, (b) Surface scanning electron 
micrograph of FeV04 thin film, (c) Cross-sectional SEM image of FeV04, thin film and (d) room 
temperature Raman spectrum on FeV04 bulk powder and thin film samples. The peak at 560 
cm~^ (indicated by an asterisk) arises from the silicon substrate 



8 



16.70 



16.65 



W 



16.60 



16.55 




10 20 30 40 50 60 70 

T(K) 




15.98 



10 12 14 16 

T(K) 



18 20 



FIG. 2: (a) Zero field temperature dependent dielectric constant (left axis) and resistivity (right 
axis) for ceramic FeV04 sample, (b) Temperature dependence of the dielectric constant near the 
magnetic ordering temperature, and (c) Temperature dependence of dielectric constant at H= 0, 
20, 40, 60, and 80 kOe. The dashed line in (b) is a guide to the eye. 



constant increases smoothly with further coohng. Since the resistivity of FeV04 increases 
monotonically with decreasing temperature (except for a small anomaly at T/V2), also shown 
in Fig. 2(a), we attribute this increase in the dielectric constant to a quartic magnetoelectric 
coupling, V4. Although FeV04does not order magnetically until cooled below T/vi=22 K, 
heat capacity measurements suggest the presence of short range spin correlations develop- 
ing well above this temperature |27j. It has been suggested in a number of other systems, 
including TeCuOs [26] and Mn304 [H], that short-range magnetic correlations can produce 
magnetodielectric corrections; we propose that the same mechanism is responsible for the 
non-monotonic temperature dependence of the dielectric constant of FeV04 in the param- 
agnetic phase. 

The fourth order magnetoelectric coupling contains terms quadratic in P and cr„(g). If 
(Ji{q) is the order parameter that develops at T/vi, then this coupling is probably domi- 
nated by A|Pp|o"i(g)p at high temperatures. The finite spin correlations developing above 
Tjvi cause (|o"i(g)p) to be nonzero so that the coupling term is in effect a|-Pp, where 
a = A(|(Ti(g)p). This term produces a shift in the dielectric constant in the paramagnetic 
phase, as seen in Fig. 2(a) below approximately 30 K. Below T/vi, approximately 22 K, when 
cri(g) acquires a finite expectation value, a trilinear coupling term aP < ai{q) > cr2(g)* is 
allowed. This term will lead to mode mixing so that the critical mode approaching T/v2 is 
not (T2{q) but {o'2{q) +pP) where p is of order a(o"i(g)) [T^. Then the divergence in this vari- 
able as T/V2 is approached will lead to a simultaneous divergence (with a very much reduced 
amplitude) in the observed dielectric constant. This mode mixing is therefore expected to 
lead to a slight increase in the magnetodielectric shift below T/vi \18\ ■ This is seen clearly in 
Fig. 2b, where the dashed line shows the extrapolation of the magnetodielectric shift from 
T > T/vi to T < T/vi, an extrapolation that does not include any corrections arising from 
the trilinear magnetoelectric term. 

C. Magnetic Field Dependent Dielectric Measurements on FeV04 Ceramic 

To further investigate spin-charge coupling in FeV04, we plot the temperature dependent 
dielectric constant measured at different magnetic fields in Fig 2(c). We find that the 
dielectric anomaly signaling the onset of ferroelectric order shifts to lower temperatures 
with increasing magnetic field, with the reduction in transition temperature reaching 0.7 
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FIG. 3: (color on-line)Magnetic field dependence of the relative change in the dielectric constant 
for ceramic FeV04 at different temperatures with vertical offset included for clarity. 

K in a magnetic field of H=80 kOe. This result is expected, as the ferroelectric order 
producing the dielectric anomaly is associated with the incommensurate spiral transition, 
which typically show a reduction in transition temperature in applied magnetic fields. 

We conducted additional measurements of the dielectric response of the bulk sample while 
sweeping the magnetic field at fixed temperature. These results are shown in Fig. 3, plotted 
as Ae(if )/e(if = 0) versus H with data measured at different temperatures offset vertically 
for clarity. At T=17 K, which is intermediate between T/vi and T/v2, there is a small negative 
magnetocapacitance, with the dielectric constant being reduced by approximately 0.03% in 
a field of H=80 kOe. As the temperature approaches the multiferroic transition at Tn2, the 
magnetodielectric coupling shows quahtative changes. By T=15 K the magnetocapacitive 
shift is positive for small fields, with a shift in dielectric constant on the order of 0.02% at 
high magnetic fields. The magnetocapacitive response is maximal near T=14.5 K, with the 
dielectric constant being reduced by approximately 0.1% in a field of H=80 kOe. At still 
lower temperatures the magnitude of the magnetocapacitive shift becomes smaller. 

Perhaps the most dramatic feature in the isothermal magnetocapacitance curves pre- 
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sented in Fig. 3 is the presence of clear maxima, which vary as a function of temperature 
and magnetic field. These maxima appear first at small fields at T=14.5 K, then shift to 
larger fields as the temperature is reduced. We believe that these anomalies do not reflect 
the suppression of the multiferroic transition temperature in a magnetic field, as discussed 
in the context of Fig 2(b). These isothermal dielectric anomalies persist to temperatures 2 
K or 3 K below T/v2, while the maximum suppression of T/v2 was only 0.7 K over the field 
range studied, as determined from the measurements in Fig. 3. We propose that this dielec- 
tric anomaly may indicate a spin-reorientation transition in FeV04. The magnetodielectric 
coupling is expected to depend on the symmetry of the magnetically ordered state P^PB] , 
so a field induced spin reorientation crossover could potentially produce the low temperature 
dielectric anomalies observed in Fig. 3. Similar magnetic field-induced dielectric anomalies 
have been observed in other materials including Mn304 [T3], although the specific mecha- 
nisms responsible remain unclear. One possibility is that the external magnetic field serves 
to reduce the slight geometrical frustration present in FeV04 [2H] , allowing a different spin 
structure to emerge. Alternatively, the spin orientation could be a spin-flop transition as 
seen in TbMnOs.jS] We note, however, that FeV04 remains ferroelectric at high magnetic 
fields |27], so the modified spin structures would still need to transform as defined by Eq. 3. 

D. Magnetoelectric Coupling in FeV04 Thin Films 

The trilinear magnetoelectric coupling that produces multiferroic order, given in Eq. 6, 
also results in an electric field (E) dependence of the magnetic structure through the coupling 
term in the free energy AF = —P ■ E. We first confirmed that these thin film samples 
were also multiferroic, through measurements of the dielectric constant and pyrocurrent, 
illustrated in Fig. 4(a). The dielectric constant for the thin film FeV04 is slightly higher than 
that found for the ceramic sample. We attribute this discrepancy mainly to the uncertainty 
in accurately determining the geometrical factor for these thin films. The dielectric response 
for these thin film samples is approximately independent of measuring frequency and the 
loss for these films is tan 6 ~ 0.01, which may be due to the presence of pinhole defects in 
the thin film sample as seen in the SEM micrograph in Fig. 1(b). The zero-field temperature 
dependent dielectric constant, measured at /=30 kHz, is plotted in Fig. 4(a). There is a 
sharp peak near T/v2=15 K, associated with the development of ferroelectric order in these 
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FIG. 4: (a) Temperature dependence of dielectric constant for FeV04 thin films at zero field, Inset: 
Zero magnetic field polarization for FeV04 thin film measured at poling fields ^poie = ±10 MV 
m~^ and (b) Temperature dependent dielectric constant measured at E=0 and E=3.75 MV m^^ 
(background was subtracted for clarity) 
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thin film samples. We note that, unlike the measurements on bulk FeV04 shown in Fig. 
2(a), the background dielectric constant for FeV04 decreases monotonically with decreasing 
temperature. This behaviour can be associated with the much larger conductivity of the 
thin film sample, arising from the presence of the pinhole defects, which obscures the low 
temperature increase in dielectric constant observed in bulk FeV04 (Fig. 2(a)). 

We confirmed that the low temperature phase of the FeV04 thin film is ferroelectric by in- 
tegrating the pyrocurrent after poling at positive and negative fields to yield the spontaneous 
polarization. These results are shown in the inset to Fig. 4(a) and indicate a spontaneous 
polarization of 6 yuC/m^, consistent with previous measurements on polycrystalline bulk 
FeV04 [27] Measurements of the dielectric response for FeV04 thin films under applied mag- 
netic fields (not shown) yield a suppression of the multiferroic transition temperature very 
similar to that observed in bulk FeV04 (see Fig. 2(b)). 

To probe the electric field control of the multiferroic phase transition temperature, ex- 
pected from the nature of the magnetoelectric coupling, we measured the temperature de- 
pendent dielectric response in the FeV04thin film sample as a function of bias voltage. 
Focusing on thin film samples allows the application of relatively large electric fields (on the 
order of MV/m) with small applied bias voltages. We chose to probe the transition through 
dielectric measurements as the magnetic anomaly at T/v2 cannot be clearly discerned in these 
thin film samples. We plot the temperature dependent dielectric constant measured at E=0 
and E=3.75 MV/m in Fig. 4(b). With the application of an electric field, the dielectric peak 
shifts upwards in temperature, by approximately 0.25 K in a field of E=3.75 MV/m. We note 
that any sample heating, which is expected to be negligible in any case because of the low 
dissipation, would raise the sample temperature relative to the thermometer temperature, 
leading to an apparent decrease in transition temperature, rather than the increase seen in 
Fig. 4(b). This increase in transition temperature is consistent with an external electric field 
promoting the development of ferroelectric order, and is similar to what has been observed 
previously in multiferroic Ni3V208 films [19]. The relatively small increase of the ferroelec- 
tric transition under such large applied electric fields can be directly attributed to the very 
small polarization in FeV04. We confirmed that the dielectric anomaly in Fig. 4(b) can still 
be associated with the multiferroic transition, even in the presence of an electric field, by 
measuring the response under the simultaneous application of magnetic and electric fields 
(not shown). Although the dielectric peak broadens considerably, the continuing presence 
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FIG. 5: (a) Electric and magnetic field dependence of the multiferroic transition temperature 
Tjv2 and (b) Magnetic field dependence of multiferroic transition temperature together with the 
proposed magnetic field induced spin reorientation cross over. Here, NCI* indicates the proposed 
phase having a spin reoriented structure 

of a single peak under such crossed fields is strong evidence that this anomaly reflects the 
multiferroic transition in FeV04. 

E. Magnetoelectric phase diagram for FeV04 

We summarize the results of these magnetoelectric and magnetodielectric studies on 
FeV04in Figs 5(a) and 5(b). We plot the E-field and H-field dependence of the multiferroic 
transition temperature in FeV04 (T/V2) in Fig. 5(a), where CI and NCI represent the incom- 
mensurate magnetic structures below T/vi and T/v2 respectively. This transition temperature 
is monotonically suppressed in an applied magnetic field, decreasing by approximately 0.7 K 
in an applied field of H=80 kOe. The transition temperature, however, increases systemat- 
ically with increasing bias voltage, shifting upwards by 0.25 K in an electric field of roughly 
4 MV/m. As discussed in Ref. [19], the magnetic field dependence of the CI-NCI phase 
boundary is expected to follow AT/v oc iJ^/^, while the electric field dependence should be 
AT/v oc E^^^^^"'\ This ability to control the transition temperature using either magnetic or 
electric fields is a key feature for a number of proposed applications for multiferroic materi- 
als. Although the size of the transition temperature shifts in FeV04 are likely too small to 
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be of any practical use, these results, taken in conjunction with previous studies on Ni3V208 
thin films, provide important evidence that this behaviour is generic among multiferroic 
materials. 

The magnetodielectric coupling in FeV04 allows us to tentatively identify the onset of 
short-range magnetic correlations, as indicated by the increase in dielectric constant below 
T=30 K in Fig. 2(a), and also to propose the onset of a spin reorientation crossover, based 
on the field dependent dielectric anomalies in Fig. 3. Using the data from Fig. 3, we 
plot this proposed spin reorientation cross-over boundary line in Fig. 5(b), together with 
the magnetic field dependence of Tjv2 (similar to that shown in Fig. 5(a)). The high-field 
putative spin-reorientation structure is labeled as NCI*. As the two boundaries do not 
coincide, the dielectric anomalies in Fig. 3 are not likely to be associated with the T/vi 
to T/V2 magnetic transition, but may potentially be attributed to a change in magnetic 
structure. Magnetic field dependent specific heat measurements (not shown) do not show 
any additional anomalies at this proposed cross-over, suggesting there is a negligible change 
in entropy between the two spin structures. 

IV. CONCLUSIONS 

We have presented a model for the development of multiferroic order in FeV04 in the 
context of Landau theory, which we have used to develop constraints on the possible mag- 
netic structures based on symmetry considerations alone. One of the noteworthy predictions 
of this model is that the electric polarization in the multiferroic phase is able to develop 
along any direction. To further investigate the higher order magnetoelectric coupling in this 
system, we have investigated the ferroelectric and dielectric response in FeV04to applied 
magnetic and electric fields. The multiferroic phase transition temperature can be tuned 
by applying electric or magnetic fields, in line with the predicted trilinear magnetoelectric 
coupling. We find evidence for a shift in dielectric constant well above the magnetic tran- 
sition temperature T/vi, which is expected to develop from a fourth order magnetoelectric 
coupling term when short range spin correlations develop in the paramagnetic phase. The 
dielectric constant shows a small, but distinct, increase below the first magnetic order tran- 
sition, which is consistent with the contribution from a trilinear magnetoelectric coupling 
term. We find evidence for magnetic field induced dielectric anomalies in the non-collinear 
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incommensurate magnetic phase of FeV04, which we attribute to a spin reorientation transi- 
tion that does not suppress the ferroelectric structure. These studies on FeV04 demonstrate 
the rich spin-charge couphng present in many multiferroic materials, emphasize the impor- 
tance to considering higher order expansions of the magnetoelectric coupling to adequately 
explain the properties of these materials, and illustrate how dielectric spectroscopy can be 
a valuable tool for probing the magnetic structures in such systems. 
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Appendix A: Landau Theory 

As discussed in detail in Ref. [18] the Fourier transform of the distribution just below 
a continuous magnetic ordering transition is proportional to the critical eigenvector of the 
inverse susceptibility matrix. (The critical eigenvector is the one whose eigenvalue first 
approaches zero, i. e. which first becomes unstable, as the temperature is lowered through 
the ordering transition.) We introduce the inverse susceptibility as follows. The thermally 
averaged spin at the site at position r in the unit cell at R, {S{R,t)) is defined as 

{S{R,t)) ^ Tr[p5op(i?,T)], (Al) 

where Sop{R, t) is the quantum spin operator at site R + t and p is the density matrix: 

p = expi-m)/[MeM-m] , (A2) 

where H is the Hamiltonian of the system. The Fourier transform of the spin distribution 
is given by 

S{q, t) = N-' Y.0iR, r))e''-' , (A3) 

R 

where r is the actual position i? + r of the spin and A^ is the total number of unit cells in the 
system. Following Landau, we write the free energy, F as an expansion in powers of S{q, r) 
as 

P=l J2 F^,r;pA<l)SMrySMr') + 0[Sm, (A4) 
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where the matrix F is the Hermitian inverse susceptibihty matrix. Of course, we do not 
know or wish to consider the exact form of "H and we do not attempt to construct the 
inverse susceptibihty from first principles. But we can analyze how symmetry influences the 
structure of the inverse susceptibility. In what follows we assume that the wave vector q at 
which ordering occurs has been established experimentally and therefore we focus only on 
that wave vector. 

We now consider the case of FeV04 which hase six spin sites within the unit cell of the 
space group PI. The only point group symmetry element is spatial inversion about the origin 
X, so that the six sites consist of three pairs of sites fn and r„+3 = — r„, with n = 1,2,3. 
Since the group of the wave vector contains only the identity element, the standard analyses 
based on this group would indicate that an allowed spin distribution function is a basis 
function of the identity irrep and therefore that symmetry places no restriction on the form 
of the spin distribution function. However, since X is a symmetry of the system when 
all the spins are zero, the free energy of the system for a configuration with an arbitrary 
distribution of Saiq, t) is the same as that for a configuration obtained by inversion applied 
to the distribution Sa{q,r). So we consider the effect of inversion on Sa{q,T). The effect 
of X is to move a spin, without changing its orientation (because spin is a pseudo vector), 
from an initial location r, to a final location — r. This means that 

I{S{R,t)) = {S{-R,t)) , (A5) 

where, for n = 1,2, 3, 

T„ = — r„ = Tn+3 = Tn , T„+3 = — r„+3 = r„ = T^^rpa • (A6) 

It then follows that 

XS{q,r) = S{q,Tr . (AT) 

Because we have 6 spins in the unit cell each having three Cartesian spin components the 
matrix F is an 18 x 18 matrix which we write in terms of 9 x 9 submatrices (for n = 1, 2, 3 
and n = 4,5, 6, respectively) as 



A B 
Bt C 



(A8) 
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Now we consider the invariance of the free energy under spatial inversion: 

= 2 E F^,,,pym^SMrr][TSp{q,r')] + 0[S{q)' 

q,a,T,l3,T' 

= l^J2 F^,r-AAq)SMr)SMry + 0[S{q)'] 

q,a,T,l3,T' 

= 2 E F^,,,^y{q)*SMrrSp{q,^') + 0[S{q)'] 

q,a,T,l3,T' 

= 1 T. F,^^,p^^iq)*SMrrSp{q,r') + 0[S{q)']. 

q,a,T,l3,T' 



(A9) 



The next-to-last equality follows because the free energy is real. The last equality is obtained 
by interchanging the roles of the dummy variables r and r and the roles of r' and r'. 



We now compare Eq. (A4) and the last line of Eq. (A9). Since these forms have to be 



equal irrespective of the values of the 5"s, we must have that 



Fa,T;P,r'iq) = i^«,r;/3,r' (<?)* • 



(AlO) 



This equality relates (for 1 < r, r' < 3) the submatrices A and C and (for 1 < r < 3 and 
4 < r' < 6) B and B"''. As a result we see that B"!" = B*, so that B is symmetric and 
C = A*. Thus 



A B 

B* A* 



(All) 



Then the eigenvectors [^,$], written in terms of the nine component vectors \1/ and $ 
satisfy 



A* + B$ = A* , 



B** + A*$ = A$ 



(A12) 



For instance 

* = [5x.(g,l), Sy{q,l), SM^), SM^). Syiq,2), 5,(g,2), S.(g,3), S,(g,3), 5,(g,3)] 
$ = [S,{qA), Sy{q,4), 5.(g,4), 5,.(g,5), 5,(g,5), 5,(g,5), 5,(g,6), Sy{q,6), 5,(g(Ap) 



The second equation of Eq. (A12) can be written as 



B** + A#* = A$* 



(A14) 
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So if [^, $] is an eigenvector with eigenvalue A, then so is exp(2p)[$*, ^*]. In principle, 
these could be two independent degenerate eigenvectors. But if one considers the simple 
case when A = aE and B = 6E, where a, and b are scalars and E is the unit matrix, one 
sees that these two solutions are, apart from a phase factor, the same. Only for special 
values of the matrices are these two eigenvectors distinct degenerate solutions. This is an 
example of an accidental degeneracy whose existence we exclude. Therefore the condition 
that these two solutions only differ by a phase factor leads to the result that 

* = e'P^* . (A15) 

Thus the nth eigenvector is 

A„, = [e'P"^l, $„] = e^^"/2[gipn/2^*^ e-'f-/^^] , (A16) 

which we write in canonical form as 

A„ = a„(g-)[e:,0„] (A17) 

where cr„(g) = |crn(^| exp(i0„) is a complex-valued amplitude and 0„ is normalized: 

1 = Ei[0„br. (A18) 

Since the inverse susceptibility matrix is 18 dimensional, there are 18 eigenvectors, each of 
this canonical form. We identify cr„(g) as the order parameter which characterizes order 
of the nth eigenvector. As the temperature is lowered one such solution (which we label 
n = 1) becomes critical and at a lower temperature a second solution (which we label n = 2) 
becomes critical. As we shall see in a moment, the magnitudes of the associated order 
parameters cr„(q) and their relative phase are fixed by the fourth order terms in the free 



energy which we have so far not considered. Using Eq. (A7) we see that 

XA„ = X[a„(g)e:,a„(g)0„] = [a„(g)*e:, a„(g)*0„] = a„(g)*An , (A19) 

which indicates that the order parameter transforms under inversion as 

Taniq) = <yn{q)* • (A20) 

Also, under spatial translation, T^, we have that 

Tj^a^{q) = e^'^-^oM ■ (A21) 
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Note that Eq. (A20) does not imply that the nth eigenvector is invariant under inversion 



about the origin. However, as we now show, it does imply that the nth eigenvector is invariant 
about an origin which depends on the choice of phase of the nth eigenvector. (It is obvious 
that a cosine wave is only inversion invariant about one of its nodes which need not occur 
at the origin.) If X^ denotes inversion about the lattice vector R, then we have 

= T^e^«-^a„(g)* = e^''^'^"'^'^ . (A22) 

Let o"„(g) = an{q)\e^^. Then if we choose R so that q-R = x, then 

XR(T„(g) = (Tn{q) . (A23) 



So, Eq. ( A20 ) implies inversion symmetry about a point which can be chosen to be arbitrarily 



close to a lattice point for an infinite system. 

Thus the contribution to the free energy from these order parameters o"„(g) at wave vector 
(f can be written as 

"■ L 

+ Yl Cnm\crn{q)o-m{q)\'^ + Yl ( ^™m[cr„(g)cr„i(g)*]^ + d*^[an(g) Vm(g)]^ j , (A24) 

where translational invariance indicates that for an incommensurate wave vector the free 
energy is a function of |o".mP, |cr„p, cr„cri^, and cr*(Tm- In writing this free energy we have 
assumed that the wave vectors of ai and o"2 are locked to be the same, as discussed in Ref. 



The generic situation in multiferroics is that as one lowers the temperature an order 
parameter o"i first becomes nonzero and then, at a lower temperature, a second order pa- 
rameter (72 becomes nonzero. In many cases, such as Ni3V208[7J or TbMn03[16j cri and 
(72 have different nontrivial symmetry. Here all the order parameters have the symmetry 



expressed by Eqs. (A20) and (A21). (The phase 02 of the second order parameter is fixed 



relative to that, 0i, of the first order parameter by the term in di2 in Eq. (A24). 



Finally, we consider the magnetoelectric coupling, V, in the free energy which is re- 
sponsible for the appearance of ferroelectricity (for which P 7^ 0, where P is the electric 
polarization). We write 



Fm + Fe + V 



(A25) 
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where Fu is the purely magnetic free energy of Eq. (A24), Fe is the dielectric potential 
which we approximate as F^ = (l/2)x^^P^, where xe is the dielectric susceptibility (whose 
crystalline anisotropy is neglected), and to leading order in cr„ 

2 

V= J2 m(^n,m,^crr,{q)am{q)* + al„^.^an{q)*(T.^{q)]P^ , (A26) 

n,m=l 7 

where n and m label order parameter modes and 7 labels the Cartesian component of P. 
Terms linear in a^ are prohibited because they can not conserve wave vector. Terms of 
order a^ or higher can exist. [35H37] The interaction V has to be inversion invariant. Since 
IP = —P and X|cr„p = |cr„p, we see that the terms with n = m are not inversion invariant 
and hence are not allowed. Thus 

V = ^Kai(g-)a2(g-)* + a;ai(g-)V2(g)]P, . (A27) 



Using XP^ = —Pj and Eq. (A20) we see that inversion invariance implies that a^ = ir^, 



where r-y is real. Then 

V = tY,r4aiiq)a2iq)* - ai{q)*a2{q)]P^ = 2 J] r^|ai(g)(T2(g)| sin(02 - 0i)Py • (A28) 

7 7 

Note that there is no restriction on the direction of the spontaneous polarization, so that 
all components of P will be nonzero. However, if the magnetic structure is a spiral, then 
the arguments of Mostovoy[31J can be used to predict the approximate direction of P. The 



result of Eq. (A28) is quite analogous to that for Ni3V208[7] or for ThMnOaflB], in that 
it requires the two modes cri(g) = exp(i0i)|(Ti(g)| and (72 (g) = exp(i02)| 0^2(9)! to be out of 
phase with one another, in other words that (pi ^ <f)2- 



If, as stated in Ref. |28l the eigenvector is not inversion invariant as implied by Eq. (A19), 
then one would conclude that the magnetic ordering transition is not continuous. However, 
the differences between the diffraction patters of the structure of Ref. |2Hland that suggested 
here are subtle enough j52] that our suggested structure seems probably the correct one. 
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